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Abstract. We investigate driven magnetic flux lines in layered type-II superconduc- 
tors subject to various configurations of strong point or columnar pinning centers by 
means of a three-dimensional elastic line model and Metropolis Monte Carlo simu- 
lations. We characterize the resulting nonequilibrium steady states by means of the 
force-velocity / current-voltage curve, static structure factor, mean vortex radius of 
gyration, number of double-kink and half-loop excitations, and velocity / voltage noise 
spectrum. We compare the results for the above observables for randomly distributed 
point and columnar defects, and demonstrate that the three-dimensional flux line struc- 
tures and their fluctuations lead to a remarkable variety of complex phenomena in the 
steady-state transport properties of bulk superconductors. 
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1. Introduction 

From a theoretical and experimental point of view, vortex matter in disordered high- 
temperature superconductors has attracted considerable attention during the past 
decades P . The possibility of practical applications of superconductivity depends on the 
maximum current density which superconductors can carry. In type-II superconductors, 
this is directly related to the flux pinning of quantized magnetic flux lines. Different 
experimental techniques such as magnetic decoration [2] , scanning Hall probes [3] , small 
angle neutron scattering jl], scanning tunneling microscope imaging [5], and Lorentz 
microscopy [6] have been utilized to capture the static structure of magnetic flux lines 
pinned by quenched disorder. These techniques provide images of the static structure of 
flux lines on the sample's top layer but do not allow mapping out the internal structure of 
the moving flux lines in response to the external current. By using neutron diffraction to 
image the flux lattice [7j, one can study the motion of these flux lines, e.g., by compiling 
a sequence of these images of vortex positions into a movie. 
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Aside from technological applications of superconductors, the statics and dynamics 
of vortices in type-II materials in the presence of quenched disorder and external driving 
force are of fundamental interest also, and have been studied quite extensively, both 
experimentally and theoretically. The presence of a small fraction of pinning centers 
tends to destroy the long-range positional order of the Abrikosov flux lattice and results 
in different spatial structures depending on the nature of the material defects. In systems 
with randomly distributed weak point pinning centers, the vortex lattice deforms and 
transforms into a Bragg glass with quasi long-range positional order at low fields [8]- 
[T3] , A vortex glass characterized by complete loss of translational order is observed if 
the pinning strength or magnetic fields are higher [33]- [18]. If correlated defects such 
as parallel columnar pins are introduced in the system, the effective pinning force adds 
coherently which results in a distinct strongly pinned Bose glass phase of localized flux 
lines characterized by diverging tilt modulus [E]-[2T]. (The term Bose glass stems from 
a mapping of the statistical mechanics of directed lines to bosonic quantum particles 
propagating in imaginary time [22], see also Ref. [23].) This type of artificial defects 
can be produced by energetic heavy ion radiation |24j . 

For vortices driven by the Lorentz force induced by external currents, analogous 
moving glass phases have been proposed. Disorder tends to inhibit flux line motion, 
and the competition between drive, pinning, and vortex interactions leads to a variety 
of different transport characteristics. At low drive, the flux lines remain pinned to 
the defects. At sufficiently large driving current, the flux lines will unbind from the 
pinning centers and start moving. At T = 0, this depinning transition constitutes a 
sharp continuous nonequilibrium phase transition; at finite temperature, this transition 
is rounded [25j [26]. In the presence of strong randomly distributed point defects, the 
ensuing moving glass is characterized by the decay of translational long-range order, the 
presence of stationary channels of vortex motion, and highly correlated channel patterns 
along the direction transverse to the motion [2?1 [28]. For weak point disorder, instead 
a topologically ordered moving Bragg glass ensues. For intermediate pinning strengths, 
one expects a moving transverse glass with smectic ordering in the direction transverse 
to the flow. When correlated disorder is present in the system, the nonequilibrium 
stationary state is predicted to be a moving Bose glass [28j [29]. 

The principal goal of the present Monte Carlo study is to characterize in detail the 
nonequilibrium steady states of interacting vortex systems in the presence of an external 
driving force and subject to a variety of configurations of strong point and/or columnar 
pinning centers by means of the force-velocity / current-voltage curve, the emerging 
spatial arrangement and vortex structures, the corresponding static structure factors, 
the mean vortex radius of gyration, and voltage / velocity noise features. This work 
is complementary to an earlier investigation that utilized the same model system and 
numerical algorithm, but addressed the weak pinning regime, and specifically focused 
on the velocity noise spectrum [30]. We shall specifically probe the effect of different 
defect configurations on the three-dimensional internal vortex structure that cannot be 
addressed in two-dimensional simulations, and is also not easily accessible in experiment. 
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Vortex line fluctuations should be expected to be prominent at elevated temperatures 
and near the depinning threshold; randomly placed point defects should further enhance 
thermal line wandering, whereas linear pinning centers tend to increase the vortex line 
tension [I]. Our results allow direct comparison with studies by other authors who 
implement a different microscopic scheme, namely Langevin molecular dynamics, to 
model the kinetics of driven three-dimensional vortex systems [31]- [31]. We remark 
that it is crucial to explore alternative simulation approaches for systems that are 
driven away from equilibrium, in order to ascertain that the results reflect physical 
properties rather than algorithmic and modeling artifacts. In addition, at low drives 
we shall explore the flux line creep mechanism via thermally activated double- kink 
configurations, and via vortex half-loops at intermediate drives [21]. Thus we aim to 
contribute to a better fundamental understanding of electrical and transport properties 
of type-II superconductors subject to various types of disorder that might aid in further 
optimization of their desired properties. 



2. Model system and numerical simulation algorithm 



2.1. Interacting disordered elastic flux line model 

We consider a three-dimensional vortex system in the London limit, where the London 
penetration depth is much larger than the coherence length. We model the vortex 
system by means of an elastic flux line free energy described in Ref. [19], see also 
Refs. [35]-[l2]. The system is composed of N flux lines in a sample of thickness L. The 
model free energy Fn (effective coarse-grained Hamiltonian) , defined by the collection of 
trajectories of the flux lines labeled by an index j on the zth layer, i.e., two-dimensional 
vectors Vj(z), consists of three components, namely the elastic energy associated with 
the line tension, the repulsive vortex-vortex interaction potential, and a disorder-induced 
pinning potential [30] : 
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Here, the elastic line stiffness or tilt modulus is given in terms of the energy scale 
e o — (0o/47rA a t,) 2 as Z\ ~ r _2 e hi(A a b/£ a b), where O = hc/2e is the magnetic flux 
quantum. The parameters A a b, £ a b, and T~ 2 are, respectively, the in-plane London 
penetration depth, the coherence length, and the effective mass ratio Mj_/M z , which 
pertains to a layered superconductor model [13]. The expression for the elastic energy 
(CQ) is valid in the limit |drj(;s:)/dz| 2 <C T 2 . The in-plane repulsive interaction between 
flux line elements is approximated as V(r) = 26qK (r/A a b), where K denotes the 
modified Bessel function of zeroth order. This function diverges logarithmically as 
r — ► and decreases exponentially for r ^> A a b- In our simulations, these vortex 
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interactions are cut off at distance L y /2 in all directions, and the system size is 
chosen sufficiently large in order that numerical artifacts due to this cut-off length are 
minimized. We model point and columnar pins through square potential wells of radius 
b with Vd(V) = — Y^k=i 0(^o — l r — r k I) as our coarse-grained defect potential. Here, 
denotes the Heaviside step function, N& the number of effective defect elements, rjf^ 
the spatial coordinates of the fcth pinning center, and U ~ y ln[l + (feo/v^ab)] 2 is the 
interpolated vortex binding energy per unit length. Finally, we add a phenomenological 
work term due to the driving force, W = — Ylf=i Jo fh ' r j( z ) m discretized form 
the model energy (TjQ) reads for the jth vortex line reads (with z*. = kbo) 

~ L/b / \ / n 2 1 L/b N v 

rp /j e i r i z t r r i z w , 1 irf\ / \ i \\\ 

¥o=2L b +2^ ^ V {H z k)-ri(z k )\) 

k=2 k=l i=l,i/j 

L/bo L/b 

+ Y^V D {r } (z k ))-J2fL-r 3 (zk) , (2) 

k=l k=l 

where A^v is the number of flux lines within the radius of the cut-off length. 

For our simulation, we chose parameter values corresponding to typical material 
parameters for YBCO (as listed in appendix D of Ref. |21j). Throughout this paper, 
the simulation lengths and energies are reported in units of the effective defect radius bo 
and interaction energy scale e (in cgs units). We have set the temperature to T = 10 
K, and used a pinning center radius and layer spacing 6 = maxjco, v^2£at>} — c o = 35 
A, anisotropy = 1/5, average spacing between defects d = 315 A= 9.0 b , and 
since the temperature is very low, an in-plane penetration depth A a b ~ Ao = 1190 
A= 34 bo, and superconducting coherence length £ a b ~ £o = 10. 5A^ 0.3 bo- Then 
e o — (0o/47rA a b) 2 ~ 1.9 x 10 -6 (in cgs units, with dimension energy /length). The 
energy scale in the first term of (j2J) is therefore £\ ~ 0.18 eo, and the pinning strength 
in the last term is Uq ~ 0.7809 €q. The free energy Fi is measured in units of eo^o- 



2.2. Monte Carlo simulation algorithm 

We employ the standard Metropolis Monte Carlo simulation algorithm in three 
dimensions with the above discretized anisotropic model free energy (j2J) [301 HO] • The 
simulations were performed in a system of size [L x , L y , L] = [-^= x 10A a b, 10A a b, 206o] with 
fully periodic boundary conditions at temperature T = 10 K or T = 0.002 e b . We chose 
the system's aspect ratio to be L x / L y = 2/y/3 to accommodate an even square number 
of vortices arranged in a triangular lattice [3U] ■ We have tested that with a penetration 
length of A a b = 35 b , the sharp cut-off interaction range of 5A a b in this system size had 
no effect on the equilibrium vortex configurations. In the absence of a driving force 
and any pinning sites, square numbers of randomly placed vortices were observed to 
arrange into a six-fold or hexagonal lattice after the system equilibrated. The average 
spacing of 9&o between defects gives a total number of 1710 lines of columnar defects. 
Each columnar defect contains 20 point defect elements, and gives a total of 1710 x 
20 point defect elements in the system. In order to be able to observe any effects due 
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number of vortices 


magnetic field B 


vortex lattice spacing 


16 


0.019 T 


101.0 b 


36 


0.043 T 


67 Ab 


64 


0.076 T 


50.5 b 


100 


0.119 T 


AO Ab 



Table 1. Magnetic fields and equivalent vortex lattice spacing for the various systems 
with different flux densities used in the simulations. The average distance between the 
1710 defects in each layer is 9.0 b . 

to these pinning centers, the maximal displacement for each random move was set to 
A = 0.25 bo- This is to guarantee that the flux lines will not move too fast and skip past 
the pinning sites. We studied the effect of maximal displacement A on the dynamics of 
vortices and found that if A was too small, the system would be trapped in metastable 
states and a much longer simulation time was required to reach equilibrium. We note 
that the limitation of flux line segment motion induces an artifact in our simulation [30J. 
This can be seen in the force-velocity curves in the form of the saturation of the curve 
at high driving force. For each move in the simulation, a displacement for the next step 
is randomly chosen from the interval between —0.25 and 0.25. The acceptance rate for 
each movement depends on the driving force: the larger the driving force, the larger the 
displacement that each step can take. However, the maximal displacement is reached 
when the driving force is increased to a finite value, which results in the saturation of 
the velocity as function of drive. We restricted our explorations to the critical regime 
where the force-velocity curve is not yet saturated. 

In this work, we randomly placed the vortex lines in the system at an initial 
high temperature T = 100 K and let them equilibrate for 50 000 Monte Carlo steps 
(MCS) in the absence of any driving force. After thus annealing the flux line system, 
the temperature was suddenly quenched to a much lower temperature T = 10 K. 
Subsequently, the external drive was applied for another 100 000 MCS to reach a steady 
state. Quantities of interest were then collected every 30 MCS for the next 250 000 
MCS, and averaged over a number of vortex lines and defect distributions. For all types 
of defect distributions, we tested that flux lines with this annealing method reached 
pinned configurations. Indeed, annealing yielded the highest critical currents in our 
simulations. In order to facilitate comparison with experimental results, we provide in 
table [T] the magnetic field corresponding to each vortex number (16, 36, 64, and 100) 
used in the simulations, as well as the lattice constants for the triangular arrays that 
would result in the absence of disorder. Velocities will be measured and listed in units 
of bo I MCS, and the driving force (current) will be given in units of eo/&o- 
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3. Quantities of interest 

3.1. Mean velocity 

The force-velocity or current-voltage curve provides important information on the non- 
equilibrium steady state of the driven vortex system. In experiments, flux lines are 
subjected to a Lorentz force in the direction transverse to the direction of an external 
current. In the simulation, flux lines tend to move in the direction of the driving force. 
The mean velocity of flux creep or flux flow is determined by the average of the total 
displacement of each flux line center of mass over a certain number of Monte Carlo steps 
(MCS): 

1 ( r cm,i(t + T) -r cmti (t)) 

V - = N-X r ' (3) 

i=l 

where N v is the number of vortex lines, r cmj ; denotes the center of mass of the zth 
vortex line, and r is set to 30 MCS. Experimentally, the average velocity of the moving 
vortices is directly related to the voltage drop across the sample from E = B x v/c. 
The driving force arises from the Lorentz force acting on the vortices when there is an 
external current applied to the sample, f = J x <p B / B. After obtaining the velocity 
for given drive, we can thus construct the force-velocity or current-voltage (I-V) curve 
and then determine an estimate for the critical driving force in each system. 

3.2. Radius of gyration 

The investigation of the three-dimensional structures in nonequilibrium steady states 
of moving flux lines is one of the main goals of this work. Hence, we are interested in 
a quantity which reflects the thermal spatial fluctuations of the flux lines. The mean- 
square displacement of flux lines measures the total displacement of the vortex center 
of mass from the beginning of the simulation. The radius of gyration is the root mean- 
square displacement of flux line elements from their center of mass. We use this quantity 
to directly investigate the effect of the different defect types on the flux line shapes; in 
discretized form, its explicit expression in the direction along the drive is 

/ 1 v-^ \ 1 I 2 

x& = [j^zz2^ Xi ^- Xc ^j ' ( 4 ) 

i,z 

with L denoting the number of layers (here, L = 20). Similarly, we can define the 
mean radius of gyration y g in the direction transverse to the drive, and perpendicular 
to the magnetic field. The total radius of gyration is r g = ^/x'i + y^. For systems 
with columnar defects, we use the maximal displacement in each flux line to determine 
the number of half-loop or double- kink excitations, see figure [1] [21] . If the maximal 
displacement in each flux line is approximately equal or greater than the average distance 
between defects, this line is considered to have a double- kink excitation. In contrast, 
for a maximal displacement greater than the size of defect but less than the average 
distance between defects, we consider this to be a half-loop excitation. 
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Super-kink 

Double-kink 

Half-loop 



Figure 1. Thermally activated half-loop, double-kink, and super-kink excitations. 



3.3. Static structure factor 

The static structure factor, which is related to the scattering cross section, reveals spatial 
symmetries and ordered structures of materials in and out of equilibrium. It is basically 
the Fourier transform of the vortex density- density correlation function (here, in d = 2 
dimensions) 

I{ q ) = J e'^-r,) ( n(ri) n{x ^ d d ri d d r2 = ^ n ^_ q y ; (5) 

where 

n(q) = J e - iq r n{r) d d r (6) 

is the Fourier transform of the local vortex density n(r). An ordered or quasi-ordered 
vortex structure such as the Abrikosov vortex lattice or Bragg glass would yield a 
periodic pattern of sharp peaks in the static structure factor. On the other hand, a 
disordered structure such as the vortex glass or vortex liquid would result in only a 
single diffuse peak at q = 0. 



3.4- Voltage noise spectrum 

The effect of disorder on the dynamics of driven vortices can also been studied by means 
of the velocity or voltage noise power spectrum which is defined by 

2 

(7) 

and directly reflects periodicity of the average velocity in a moving vortex lattice. The 
voltage noise spectrum can display broadband or narrowband noise. It has been shown 
that in the presence of weak point defects, and at low driving force, vortices are in 
a plastic flow regime, characterized by a broadband noise signal S(uj) ~ u~ a [44]. 
This power law may be interpreted as a remnant of the zero-temperature continuous 
depinning transition. 

Narrowband noise is expected to exist provided the moving vortices maintain a long- 
range or quasi long-range positional order as in the moving Abrikosov vortex lattice or 
moving Bragg glass. For weak point or columnar pinning centers, a thorough study of 
the voltage noise power spectrum was reported in Ref. ]30j. There, narrowband noise 



v(t)e lult dt 
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emerged at the characteristic washboard frequency which is given by u = 2ir(v)/a, with 
the average velocity (v), and the typical intervortex distance a. The reason for the 
existence of this washboard noise is that quenched disorder in the system occasionally 
traps the vortices. As a consequence of the periodic vortex structure, this results in a 
periodically varying average overall velocity of the moving Abrikosov lattice or Bragg 
glass. 

The existence of broadband noise and appearance of the washboard frequency 
in narrowband noise has been experimentally demonstrated [I5]-[IE] and numerically 
investigated both in two-dimensional [19]- [51] and three-dimensional systems [30, 33J. 
The magnitude of the broadband noise strongly depends on the applied current and the 
external field [52]. It was shown that the broadband noise in the 'peak effect' regime 
appears with the onset of vortex motion. The noise power increases to its maximum 
slightly above the critical current, and then decreases at large applied currents where the 
noise becomes suppressed. It was also suggested that the magnitude of the broadband 
noise would reach its maximum in the regime of plastic flow, and in three-dimensional 
simulations, narrowband noise with washboard frequency peaks was reported in the 
moving Bragg glass in the presence of weak pinning disorder [30]. In a recent experiment, 
narrowband noise was detected even in the peak effect regime, indicating apparent 
long-range temporal correlations of vortices near the critical temperature T c [53]. The 
characteristic frequency found in this experiment did not coincide with the washboard 
frequency, however: the associated length scale matched the sample width rather than 
the intervortex distance. Voltage noise has also been used to study the transition from 
a moving disordered phase, such as the moving vortex glass, to a moving ordered phase, 
such as moving Bragg glass or moving Abrikosov vortex lattice [33, 49J. The existence 
of washboard noise in systems with strong point or columnar defects is one of the issues 
to be addressed in this present study. 

4. Simulation results 

4-1- I-V characteristics 

We first address the effect of different defect configuration on the I-V curves. In 
our simulations, we have studied six defect configurations: regular triangular and 
rectangular arrays of columnar defects; randomly placed splayed columnar defects; a 
mixture of point and randomly placed columnar defects; and randomly placed point 
defects. In order to generate columnar pinning centers, defect elements in each layer 
are arranged on top of each other which results in highly correlated defects with the 
highest critical currents. For our system with mixed point and linear defects, we chose 
the ratio between point and randomly placed columnar defect elements as 1:1. Splayed 
defects are created by placing columnar defects at random positions and tilting them 
slightly at fixed angle in random directions. In our simulations we tilt them such that 
the displacement between the top and the bottom layer for each defect is about 5 bo. 
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Figure 2. Small sample cross sections of systems with (a) rectangularly, (b) 
hexagonally, and (c) randomly placed columnar defects, (d) a mixed system with both 
randomly placed point and columnar pinning centers, and systems with (e) randomly 
placed splayed linear pins, and (f) point defects. 



(A larger lateral displacement value would result essentially in a randomly placed point 
defect configuration and therefore lower critical currents.) With the same defect density, 
the difference in the value of critical currents for each I-V curve directly reflects the 
distinct spatial defect configurations. Small sample slices for each defect configuration 
type are depicted in figure [2j 

Characteristic I-V curves for systems with different types of defect distributions 
are shown in figure [3j It is apparent that the critical currents for systems with 
triangular and rectangular arrangement of columnar defects are considerably higher 
than in samples with randomly placed columnar pins, randomly placed splayed columnar 
defects, a mixture of point and randomly placed columnar pinning centers, and randomly 
placed point defects, respectively. This is in good agreement with the experimental and 
numerical work [TJ [2U [30] , which clearly demonstrate that columnar defects are more 
effective and thus yield a larger critical current density J c than point pins. At low 
and intermediate drive, flux lines are trapped in metastable states and spend most 
of their time inside or in the immediate vicinity of the pinning potential. Flux line 
wandering between nearby defects may occur due to the thermally activated jumps. 
The average velocity is extremely small and can be explained by the theory of flux creep 
[SI]. Far above the critical current, the flux lines flow freely. In the case of a system with 
randomly placed point defects, the pinning force density does not add coherently over the 
length of the vortex. Instead, the flux lines will attempt to find energetically optimized 
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Figure 3. Force- velocity or FV curves for systems with different defect configurations. 
Each system contains 16 flux lines. Fhe systems with rectangularly (circles) and 
triangularly (squares) arranged columnar defects have higher critical currents than 
samples with randomly placed columnar defects (diamonds), splayed columnar defects 
(left triangles), a mixture between points and columnar defects (up triangle), or point 
defects (down triangle), respectively. This confirms that systems with correlated 
defects yield higher critical currents. Note that here the pinning strength (per unit 
length) for different defect configuration is set to be the same, 0.78 eobg, while the real 
pinning strength for point defects is smaller by approximately an order of magnitude. 
Each of the data points was obtained from taking an average over 50 different defect 
realizations. Driving forces are given in units of and vortex velocities are 

measured in units bo / MCS. 

paths through the point pinning landscape, which promotes flux line wandering in the 
sample. Note that the real pinning strength for a columnar defect is stronger than for 
a point defect by approximately an order of magnitude. In our simulations, the pinning 
strength (per unit length) of all defect configurations is set to be the same, comparable 
to typical columnar pin strengths, whereas in Ref. [30] a much weaker pinning potential 
was chosen, more in line with typical oxygen vacancy pinning strengths. We have also 
confirmed that the depinning current (critical current at T = 0) is lowered for all types 
of defects in denser vortex systems, where the repulsive vortex interactions tend to 
dominate over the attractive pinning energies [301 I3T1 [32] . This is demonstrated in 
figuress H] and [5] for randomly placed columnar and point defects, respectively, showing 
the I-V characteristics obtained for systems with 16, 36, 64, and 100 lines in either case. 
In figure El we compare the approximate depinning forces or critical currents for systems 
with varying flux densities and randomly arranged columnar or point pinning centers, 
as obtained from the data shown in figures H] and El 

4-2. Radius of gyration 

Next we determined the root mean-square displacement or radius of gyration PJ) near 
the depinning current, which allows us to obtain additional information about the shape 
of the moving flux lines in our three-dimensional samples. (This quantity is obviously 
not available in a two-dimensional simulation.) The radius of gyration represents the 
spatial fluctuations of the flux line from its center of mass. We omitted results in the 
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Figure 4. Force-velocity or FV curves for systems with 1710 randomly placed 
columnar defects and different flux densities: 16 lines (circle), 36 lines (square), 64 
lines (diamond), and 100 lines (triangle up). Fhe depinning force or critical current 
becomes smaller for denser vortex systems. 




Figure 5. Force-velocity or I-V curves for systems with 34200 randomly placed point 
pinning centers and different flux densities: 16 lines (circle), 36 lines (square), 64 lines 
(diamond), and 100 lines (triangle up). Again, the depinning force or critical current 
goes down with increasing vortex density. 



(a) 
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Figure 6. Critical current or depinning forces for (a) systems with randomly arranged 
columnar pins and (b) point defects as function of vortex density, as obtained from the 
data shown in figures [4] and [5] 
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Figure 7. Mean radius of gyration (in units of bo) along the drive direction for systems 
with randomly placed columnar defects with different flux densities: 16 lines (circle), 36 
lines (square), 64 lines (diamond), and 100 lines (triangle up). The radius of gyration 
tends to grow at all drives as the density of flux lines increases, but saturates for dense 
systems. Recall that the mean columnar defect separation is 9 bo- 



a 1 s 

2 1.2 
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Figure 8. Mean radius of gyration in the direction transverse to the drive for systems 
with randomly placed columnar defects for different flux densities: 16 lines (circle), 36 
lines (square), 64 lines (diamond), and 100 lines (triangle up). 



regime of low driving forces since the dynamics there is extremely slow. Data in the 
regime of large driving forces are discarded also due to the unphysical limitation for 
each line element displacement as mentioned in the previous section. 

4-2.1. Columnar defects. Figures [7] and El display our results for the average radius of 
gyration (Tj0) as functions of the driving force in systems of 16, 36, 64, and 100 flux lines 
with 1710 randomly placed columnar defects. For all investigated vortex densities, the 
radius of gyration in the direction of the drive increases with the driving force in the 
range studied here. At low drives, we observe the Bose glass phase with completely 
localized vortex lines inside the columnar defects as shown in the snapshot in figure [9h- 
In this state, the mean radius of gyration would be smaller than the columnar pin radius, 
which is set to 1 in our simulations. Since the pinning force of a columnar defect is very 
large and overcomes the repulsive vortex interaction, the flux lines become localized at 
different defect locations resulting in the disordered Bose glass phase. At intermediate 
drives close to the critical regime, some flux lines become delocalized and move along the 
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(a) (b) H§ 



Figure 9. Snapshots of (moving) vortices in the presence of randomly placed columnar 
defects at (a) low, (b) intermediate, and (c) large drive. The flux line system is in the 
(a) pinned Bose glass phase, (b) partially pinned Bose glass and moving vortex glass 
/ liquid phase, and (c) moving vortex glass / liquid phase. 

drive direction, with enhanced transverse line fluctuations. This coexistence of a fluid of 
moving flux lines and the localized Bose glass is shown in figure [Ub. One would expect 
a small average radius of gyration induced by the moving flux lines. At high drives, 
all flux lines are in motion, and hence a larger radius of gyration ensues, see figure [9b. 
The radius of gyration also increases with the number of vortices, but saturates at large 
flux densities, in our system at 64 and 100 lines, as shown in figure [7J This increase 
with flux density is again a consequence of the repulsive vortex interactions beginning 
to dominate over the attractive pinning energies. Flux lines are increasingly caged 
by their neighbors, and transverse wandering away from the linear defects becomes a 
collective excitation affecting several vortices. The observed saturation indicates that 
caging due to the effectively stronger vortex interaction at large densities governs the 
flux line dynamics at high driving currents. Notice that for the largest flux density in 
our simulation, the gyration radius reaches beyond the mean columnar pin distance, 
allowing for pinning of single lines to several defects. Qualitatively similar behavior of 
the mean radius of gyration is seen both along and transverse to the drive direction; 
vortex line fluctuations along the drive direction are however an order magnitude larger 
than the transverse fluctuations. 

Note that the above results differ characteristically from the data reported in 
Ref. [30] for systems with much smaller pinning strengths. In the moving vortex states 
of the weak pinning regime, the radius of gyration, both along the drive and transverse 
to it, was found to decrease for increasing flux density. This situation corresponds 
to a dynamics which is dominated by the vortex interaction rather than the pinning 
energies. This marked difference to the trend reported for our presented data arises 
from the fact that the pinning strength used in the current simulations is stronger by a 
factor of approximately 10: here we address a defect-dominated regime. In addition, we 
simulated our system with a larger penetration length which requires a larger system 
size in order to avoid the numerical artifact due to the interaction cut-off. With the same 
number of vortex lines, but in a larger system size and much stronger pinning strength, 
we are here working in a regime where the flux lines are rather weakly interacting, 
and their structure and dynamics is dominated by the defects. Thus, the investigated 
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Figure 10. Mean radius of gyration along the drive direction for systems with 
point defects with different flux densities: f6 lines (circle), 36 lines (square), 64 lines 
(diamond), and 100 lines (triangle up). The radius of gyration starts to grow as the 
drive increases, reaching its highest value at the critical current; at larger drives, it 
decreases again. Note that even for the largest flux density the peak gyration radius 
is only about 1/2 the corresponding vortex lattice constant, see table [T] 
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Figure 11. Mean radius of gyration transverse to the drive and magnetic field for 
systems with point defects with different flux densities: 16 lines (circle), 36 lines 
(square), 64 lines (diamond), and 100 lines (triangle up). The radius of gyration 
grows monotonously as the drive increases. 

current range effectively remains rather close to the depinning threshold. 

4-2.2. Point defects. Point defects directly affect the vortex dynamics. In contrast to 
systems with columnar defects, uncorrelated point disorder promotes spatial wandering, 
transverse to the magnetic field direction, of the moving flux lines. The spatially 
randomly distributed pinning centers try to pull the vortex lines in different directions 
as they traverse the sample, which leads to less effective pinning compared to linear 
defects. In the regime where the defects dominate the dynamics, the radius of gyration 
for point defect is then expected to be larger than for systems with columnar pins. 

Figure [10] shows the radius of gyration along the direction of the drive. The shape 
of the curve differs from the one depicted in figure [7] for randomly placed columnar 
defects. At small drives, the radius of gyration at a specific driving force slightly 
increases with growing flux density, while we observed a larger deviation in the system 
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Figure 12. Snapshots of moving vortices in the presence of randomly placed point 
defects at (a) low, (b) intermediate, and (c) large drive. Flux lines are in the (a) pinned 
vortex glass phase, (b) and (c) moving liquid / vortex glass phase. The vortices display 
the largest average radius of gyration at intermediate drive. 

with columnar defects. For increasing drives, the radius of gyration keeps growing and 
reaches its maximal value just above the depinning threshold. This indicates that the 
dynamics of the system is only marginally dominated by the point pinning centers, and 
only until the critical threshold is reached. In an infinite system at zero temperature, 
critical fluctuations would lead to a divergence of the mean radius of gyration at the 
nonequilibrium depinning phase transition. In the moving glass phase, we observe that 
the gyration radius decreases with flux density. At lower densities, point disorder plays 
a more prominent role, and has the effect of enhancing flux wandering which allows the 
vortex lines to more efficiently explore favorable sites in the pinning landscape. We note 
that the gyration radius along the drive direction peaks at values well below the typical 
vortex lattice constants, see table [H 

The radius of gyration in the direction transverse to the drive behaves quite 
differently. As shown in figure [Tl] at any vortex density the transverse radius of gyration 
keeps monotonously increasing with the driving force, and no tendency of a maximum 
or saturation is evident. The curves tend to merge towards the same value at high 
drive. The transverse radius of gyration is hardly affected by the drive; flux lines can 
only diffuse about their center of mass in the absence of the drive. The distance that 
the flux line element can move away from the center of mass is limited by the elastic 
line tension, which tries to pull the vortex line element back towards its center of mass. 
For increasing vortex density, the radius of gyration becomes larger since the dynamics 
becomes increasingly dominated by the vortex interactions rather than the pinning 
potentials; displacement of a line element of any vortex thus typically induces motion of 
nearby elements of other flux lines as well. Snapshots of moving vortices subject to point 
pinning centers at different drives are plotted in figure [121 Figure [12b at intermediate 
driving force corresponds to the largest radius of gyration observed in our simulations. 
Indeed, transverse line wandering is prominent near the depinning current, as is evident 
in figure [12b. (Note, however, that the layer-to-layer displacements still range within 
the applicability limit of the London approximation.) 
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Figure 13. Measured number of half- loops for systems with random columnar defects 
with different flux densities: 16 lines (circle), 36 lines (square), 64 lines (diamond), 
and 100 lines (triangle up). The number of half- loops shows a maximum just above 
the critical current and then tends to decrease at all drives as the density of flux lines 
increases. 
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Figure 14. Measured number of double-kinks for systems with random columnar 
defects with different flux densities: 16 lines (circle), 36 lines (square), 64 lines 
(diamond), and 100 lines (triangle up). The number of double-kinks tends to increase 
with growing vortex density at all drives, similar to the mean radius of gyration. 



4-3. Half-loop and double-kink excitations 

The crossover of the dynamics from a disorder- to an interaction-dominated regime 
as induced by the driving force has been explained in the previous section. A large 
mean vortex radius of gyration appears at elevated driving forces for systems with both 
randomly placed point and columnar defects. At low driving force, the structure of the 
flux line system is different in both cases. Columnar defects tend to straighten the flux 
lines, whereas in a system with point defects flux line wandering is promoted. At low 
driving force, we observed the localized Bose glass in systems with columnar defects. At 
low temperatures and low external driving current, flux creep happens via the formation 
of vortex half- loop, double- kink, and superkink excitations (see figured]) pi| [55} 156] . 
which represent saddle point configurations whose excitation energies U^{J) represents 
the effective current-depending barrier energy that enters the thermal activation factor 
exp[— [/b( J) /T] for flux flow. In the localized Bose glass phase, these excitations occurred 
extremely rarely in our simulations, preventing us from obtaining statistically significant 
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Figure 15. Static structure factor plot of moving vortex lines in a system with 
randomly placed point defects in the nonequilibrium steady state. The strong pinning 
centers destroy any translational order, and the ensuing single structure factor peak 
indicates a moving plastic or liquid phases of vortices, as also seen in the snapshot 1 12b . 
We obtain similar static structure factors at all vortex densities, driving forces, and 
for all defect configurations studied here. 

data. 

At intermediate driving forces near the depinning threshold, however, half-loop 
excitations become quite prominent, as seen in figure [131 an d at given vortex density, 
their typical number decreases with increasing drive. Also, a larger number of half-loop 
excitations is observed as the density is increased, which indicates again that the vortex 
interactions play an important role for the existence of multiple, 'coherent' half-loop 
excitations. 

We do not observe double kinks at small driving currents; they only appear at drives 
larger than the critical depinning force. The presence of double-kinks is largely due to 
an interplay of pulling on the vortices and strong pinning forces rather than thermal 
fluctuations. As a flux line is moving, some vortex elements will enter columnar pinning 
centers. These flux line elements become trapped while other segments are pulled along 
the drive direction. The vortex line keeps stretching until the elastic energy overcomes 
the pinning energy, whereupon the trapped flux line elements depin from the defect. This 
would result in the emergence of double-kinks, as depicted in figure [TH These prominent 
double-kinks are likely largely responsible for the observed continuing increase of the 
mean radius of gyration, see figure [71 

4-4- Static structure factor 

In this present study, the abundant point or linear pinning centers are always sufficiently 
strong to destroy any spatially ordered structures in the accessible range of driving 
forces. In contrast with our earlier investigations of vortex transport in the presence 
of much weaker pinning sites [30J, here we observe neither positional nor orientational 
long-range order. In this strong pinning situation, flux transport happens via plastic or 
liquid motion at drives larger than the depinning force. Correspondingly, we observe 
a single peak in the static structure factor of the moving vortex system, for all defect 



Nonequilibrium steady states of driven maanetic flux lines 



18 




to 



Figure 16. Voltage noise along the drive direction for systems with randomly placed 
point defects at various vortex densities. The driving force along the x-axis for each 
curve was set to 1.5. The noise amplitude grows as the vortex density increases from 
16 to 36, 64, and 100 lines in the system. The broadband noise in these systems reflects 
the presence of a moving disordered phase, such as a moving liquid or plastic phase. 
This is also supported by the structure factor plot and the snapshot shown in the inset. 
The power laws S(w) oc and S(ui) oc 1/u 1 - 7 are indicated also, to match to the 

data at low and intermediate frequencies, respectively. 

configurations investigated here. This corresponds to a single peak in the diffraction 
plot as the experimental signature for plastic motion j32j EZj- As an example, figure fT5l 
shows the structure factor ([5]) obtained from our simulations for a system with point 
defects. The single peak indicates the presence of an isotropic liquid or a disordered solid 
of moving flux lines, as indeed confirmed by the snapshot depicted in figure [12b . (We 
have also measured the static structure factor in a different system with weak pinning 
defects, and instead observed hexagonal Bragg peaks corresponding to the triangular 
moving vortex lattice or Bragg glass.) 

4-5. Voltage noise spectrum 

Finally, we report and discuss the measured voltage noise power spectrum ([7]) stemming 
from the vortex motion in systems with randomly placed point and columnar defects. 

4-5.1. Randomly placed point defects. As shown in figure fT6l we observe broadband 
noise in our simulations of systems with point defects at all vortex densities. The results 
are averaged over 50 defect realizations and initial vortex positions. We specifically 
investigate the regime where the driving force is set to 1.5, and after the system has 
reached its nonequilibrium steady state. The absence of narrowband noise indicates 
the presence of a moving vortex liquid phase, in agreement with experiments in the 
peak effect regime [52]. This interpretation is supported by the appearance of the single 
peak in the structure factor and the snapshot of the moving vortices in the inset of 
figure dEl Each flux line in this system subject to strong point disorder moves largely 
independently, which implies the loss of any washboard signal in contrast with the results 
obtained in samples with very weak point defects [31)] • In the presence of weak pinning 
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Figure 17. Voltage noise in the transverse direction to the driving force (y-direction) 
for systems with randomly placed point defects at various vortex densities. The driving 
force along the x-axis for each curve is set to 1.5. Similar to the voltage noise in the 
x-direction, the magnitude for each curve decreases as the vortex density increases 
from 16 to 36, 64, and 100. The broadband noises in these systems reflect the property 
of a moving disordered phase. 



centers, the vortices arrange themselves into the quasi long-range ordered moving Bragg 
glass, and essentially move with the same speed. As this ordered structure traverses the 
disordered substrate with almost homogeneous speed (v), periodic trapping induces a 
narrowband washboard noise signal with a characteristic peak frequency u = 27r(v)/ao, 
where ao denotes the vortex lattice constant. 

The magnitude of the broadband noise signal decreases for increasing vortex density, 
indicating that the flux lines tend to reorder at higher flux density. This is also reflected 
by the decrease of the radius of gyration for increasing vortex density. The magnitude 
of the noise indicates the fluctuations of the velocity about its mean. At higher flux 
densities, the vortex interactions dominate and cause a stiffening of the flux lines as 
they become caged by their neighbors. Stronger vortex interactions favor the Abrikosov 
vortex lattice or Bragg glass with reduced noise. A decrease of the noise power for 
increasing vortex density or interaction strength has also been reported in experimental 
data [52] and numerical simulations [19] . 

A power law broadband noise spectrum S(u) ~ l/u a was also reported in 
simulations of two- and three-dimensional XY and dual Coulomb gas models [H] , where 
at zero magnetic field an exponent a = 1.5 was found, whereas a ~ 1 in a system with 
small applied magnetic field. These results are similar to the experimental data reported 
in Refs. [58j [59] , where the voltage noise spectrum decays at high frequency like 1/u 
at elevated temperature and like at lower temperature. The exponent a = 1.5 

also emerges from a functional renormalization group calculation for point defects to 
one-loop order in three dimensions [60]. The voltage noise power spectrum in our data, 
with strong pinning centers, varies approximately like S(u) ~ l/u> a with a ~ 1.3, which 
is well fitted in the low frequency regime, whereas we find a ~ 1.7 at intermediate 
frequencies. Our exponent value is smaller than the exponent a = 2, which is in fact 
the mean-field value, and was measured in a system of non-interacting driven vortices 
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Figure 18. Voltage noise power spectrum in the direction of the driving force 
for systems with randomly placed columnar defects at various flux densities. The 
driving force along the x axis for each curve is set to 2.7. The amplitude for each 
curve increases with the vortex density. The low-frequency oscillations originate from 
repeated traversal through the finite system with periodic boundary conditions. 

subject to weak point and columnar defects [40] . 

As shown in figure \TT\ the y-component of the broadband noise is observed to 
be nearly flat white noise, and to decrease in magnitude for increasing vortex density. 
Transverse to the drive direction, there are much smaller velocity fluctuations, and no 
typical time scale exists. Each moving flux line can only fluctuate around its center of 
mass which results in zero average velocity. 

4-5.2. Randomly placed columnar defects. In our system with 1710 randomly placed 
parallel columnar defects, we observe narrowband voltage noise peaks on top of a 
broadband background, as shown in figure [18] for systems with a 16, 36, 64, and 100 
flux lines, with the data taken in the nonequilibrium steady state at driving force 2.7 
along the x direction, and the average taken over 50 realizations of the defect and initial 
vortex positions. As can be readily confirmed by computing the associated characteristic 
time, the presence of these narrowband peaks at low frequencies can however be simply 
attributed to the traversal time of the vortices through the entire system. The observed 
periodicity is merely due to the moving vortices encountering the identical linear defect 
configurations in our system with periodic boundary conditions. Even though our 
system with point defects has similarly strong pinning centers, no narrowband peaks 
are observed here since each flux line element moves almost independently. During their 
repeated traversal through the sample, each flux line element may become trapped at 
point pins at different locations. The various associated time scales yield quite different 
frequencies and result in the characteristic broadband noise. In comparison to the system 
with point defects, our sample with randomly placed columnar pins yields a power law 
S(u) ~ 1/oj in the low frequency regime. Thus, different defect configurations induce 
distinct values of the broadband noise scaling exponent: correlated defects produce a 
larger exponent value in the strong pinning regime. 
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5. Discussion and conclusion 

In this paper, we have presented numerical investigations of the nonequilibrium steady 
states of driven magnetic flux lines in type-II superconductors subject to various 
configurations of strong point or columnar pinning centers. We have employed a versatile 
three-dimensional Metropolis Monte Carlo code based on an elastic line model and used 
the measured current-voltage curve, vortex structure factor, mean radius of gyration, 
number of half-loop and double-kink excitations, and the voltage noise power spectrum 
to thoroughly characterize the ensuing moving vortex structures and fluctuations. 

Interesting features arise due to the competition between the elastic energy, pinning 
potential, repulsive vortex interactions, and the driving force. In the absence of the 
driving force and defects, we observed the Abrikosov vortex lattice as expected. A 
pinned vortex glass is found in systems with strong point defects, while a localized Bose 
glass is observed in the presence of strong columnar defects. From our simulation results, 
the snapshots of these two systems from the top view look very similar, i.e., each flux 
line looks like a point-like particle located inside the pinning center. However, these flux 
line systems behave quite differently if a small driving force is applied. The effect due to 
different defect structures is reflected in their different force-velocity or I-V curves. We 
have found that systems with correlated defects such as parallel or splayed columnar 
defects yield higher critical depinning currents. The presence of strong point defects 
strongly promotes flux line wandering, which results in a liquid-like vortex structure. 
This feature is also captured by measuring the radius of gyration and studying snapshots 
of the driven vortex system. By investigating the radius of gyration for various driving 
forces, we observed that the driving force changed the moving flux line dynamics from 
a regime dominated by the disorder at low drives to a regime dominated by the vortex 
interactions at high drives. This is, for example, evident in the system with point defects 
as the driving force or the vortex density are varied. In contrast, in systems with weak 
pinning centers the disorder-dominated regime is not accessible. 

In the low driving force regime, we did not observe double-kinks excitations, but 
a small number of half-loop excitations were detectable at intermediate driving forces. 
At large driving force, we observed a liquid-like or amorphous disordered structure of 
moving flux lines in our samples with strong point or columnar defects. The static 
structure factor in all our results show a single peak which corresponds to a disordered 
structure. Another noticeable difference between the vortex systems subject to point 
and columnar pinning centers emerges when studying the voltage noise spectrum. While 
just above the critical threshold broadband noise is observed at all vortex densities and 
driving forces in either case, the characteristic exponent describing the power-law decay 
of the voltage noise signal turns out to be larger for randomly distributed point defects. 
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